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Abstract 

We study gluon scattering amplitudes in N = 4 super Yang-Mills theory at 
strong coupling via the AdS/CFT correspondence. We solve numerically the dis- 
cretized Euler-Lagrange equations on the square worldsheet for the minimal surface 
with light-like boundaries in AdS spacetime. We evaluate the area of the surface for 
the 4, 6 and 8-point amplitudes using worldsheet and radial cut-off regularizations. 
Their infrared singularities in the cut-off regularization are found to agree with the 
analytical results near the cusp less than 5% at 520x520 lattice points. 



1 Introduction 



Recently Alday and Maldacena [T] proposed the prescription for the calculation of gluon 
scattering amplitudes in M = 4 super Yang-Mills (SYM) theory at strong coupling via the 
AdS/CFT correspondence. They showed that the gluon scattering amplitude is related to 
the area of a minimal surface in AdS spacetime surrounded by the Wilson loop with light- 
like boundaries. They found the exact solution of the minimal surface corresponding to 
the 4-point amplitude and showed that it reproduces the perturbative formula conjectured 
by Bern, Dixon and Smirnov (BDS) [2]. 

The correspondence between gluon amplitudes and the Wilson loops has been exam- 
ined also at weak coupling [3j. It is shown in [4] that the Wilson loops at weak coupling 
obey the anomalous conformal Ward identity, which determines the n = 4 and 5-point 
amplitudes completely. For n > 6-point amplitudes, however, the conformal invariance of 
the amplitudes does not fix their dependence on the kinematical variables. Recently it is 
found that the 6-point 2-loop corrections to the Wilson loop agrees numerically with the 
gluon amplitudes but they differ from the BDS formula[5]. 

In order to study gluon scattering amplitudes at strong coupling using the AdS/CFT 
correspondence, we need to find the solution of the minimal surface in AdS spacetime 
surrounded by the light-like Wilson loop. See [6], [TJ, [HI El El EH EEl EE] for references. The 
minimal surface is obtained by solving the Euler-Lagrange equations with the Dirichlet 
boundary conditions. These are non-linear partial differential equations and highly non- 
trivial to solve analytically for the polygon with n > 5 boundaries. Any exact solution 
for n > 5 is not yet known so far. In a previous paper [i9], the authors including the 
two authors in the present paper, constructed solutions for the 6 and 8-point amplitudes 
by cutting and gluing the 4-point amplitude. The evaluation of the amplitudes shows 
that they do not agree with the BDS formula both for the infrared singularity and finite 
parts. In particular, the infrared singularities coincide with those of the 4-point ampli- 
tudes, which suggests that these simple solutions are not the minimal surfaces and might 
correspond to the other disconnected amplitudes. 

In this paper we will propose a practical approach to compute the minimal surface in 
AdS spacetime. We investigate numerically the solution of the Euler-Lagrange equations 
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for the minimal surface in AdS spacetime, surrounded by the light-like segments. We solve 
the discretized Euler-Lagrange equations on the square lattice with the Dirichlet boundary 
conditions. We evaluate the area of the surface for the 4, 6 and 8-point amplitudes using 
two types of cut-off regularizations. One is the world-sheet cut-off regularization. The 
other is the radial cut-off regularization [13]. For the 4-point amplitude we see that the 
numerical result agrees with the analytical solution obtained by Alday and Maldacena. For 
the 6 and 8-point solutions, we take the same momenta configuration as in [9] . The results 
are different from the cut and glue solutions and there appear new cusp singularities. Their 
infrared singularities are found to agree with analytical solutions near the cusp. This 
numerical approach would be a useful method to calculate the n-point gluon amplitudes 
at strong coupling and test the BDS conjecture from the AdS side. 

The paper is organized as follows. In Section 2 we describe the discretized Euler- 
Lagrange equations in AdS spacetime and solve the equations numerically. In Section 
3, we evaluate the area of the minimal surface for 4, 6 and 8-point amplitudes by using 
the cut-off regularization in the radial direction of AdS spacetime [13] . We compare their 
infrared singularity part with the analytical results near the cusp. Section 5 includes 
conclusions and some comments. 

2 The Euler-Lagrange equations of the minimal sur- 
face 

2.1 Analytical solutions 

In this paper we investigate the minimal surface in AdSs surrounded by the curve C n made 
of light-like segments Ay 11 = 2np^ 1 which corresponds to the n-point gluon amplitude with 
on-shell momenta pi (pf — 0, i — 1, • • • ,n). Here y^ (/i = 0,1,2,3) together with the 
radial coordinate r are the Poincare coordinates in AdS^ spacetime with the metric 

ds > = R ^ d y^ +dr \ (2 .i) 

and R is the radius of AdSs. It is convenient to write the Nambu-Goto action in the static 
gauge where we put y 3 = and parametrize the surface by y\ and y 2 . Then r and y are 
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functions of 7/1 and y 2 and the action is given by 




The Euler-Lagrange equations become 

*(a&)=°- *fe)-§F = °- ^ 

where L is the Lagrangian of the action. The minimal surface is obtained by solving the 
Euler-Lagrange equations (12.31) . which are non-linear partial equations and difficult to 
solve analytically. 

For n = 4, Alday and Maldacena found the exact solution pp. Let us consider the 
scattering amplitude for two incoming particles with momenta p\ and p% and outgoing 
particles with momenta P2 and p±. The Mandelstam variables are given by s = —{pi +P2) 2 
and t = — (p 2 + Pz) 2 - The case with s = t is particularly simple where the Wilson loop is 
represented by the square with corners at 2/1,7/2 = ±1- The momenta in the (y , 7/1, y 2 )- 
space associated with the square are 

2ir Pl = (2, 2, 0), 2np 2 = (-2, 0, 2), 2ixp z = (2, -2, 0), 2n Pi = (-2, 0, -2). 

(2.4) 

The boundary conditions for the case with s = t are given by 

r(±l, 7/2) = r{ Vl , ±1) = 0, 7/ (±l, 7/2) = ±2/2, 2/0(2/1, ±1) = ±Vv (2.5) 

The solution of the Euler-Lagrange equations satisfying the boundary conditions is 

2/0(2/1,2/2) = ym, r(y u y 2 ) = \J (I - - yf). (2.6) 

By applying the conformal transformation 50(2,4) one can obtain the solution with 
general s and t. Using the dimensional regularization, the area is shown to agree with the 
BDS formula at strong coupling. 

Noticing that the solution f)2.6p satisfies eq. (12.31) even when we change the sign of 
7/0, an obvious and simple generalization of this remarkable 4-point solution is to cut and 
glue the solution. In [9J, we have constructed 6-point and 8-point solutions of the same 
action. These cut and glue solutions are summarized as follows: 
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6-point function solution 1: 

Vo = ^ CI2/12/2 1 + ym - \yi\y* H- 2/1I2/2I), r = yj(l -y\){l -yl). (2.7) 
The solution corresponds to the scattering with momenta 

2vr Pl = (2,0,-2), 2tt P2 = (-1,1,0), 27rp 3 = (1, 1, 0), 

27rp 4 = (-1,0,1), 2tt P5 = (1,0,1), 2 m = (-2,-2,0). (2.8) 

6-point function solution 2: 



2/o = 2/1I2/2I, r= ^/(l-y?)(l-2/2)- (2-9) 



The momenta are 



2tt Pi = (1,0,-1), 2np 2 = (-1,0,-1), 2tt P3 = (2,2,0), 

27TP4 = (-1,0,1), 27Tp fi = (1,0,1), 27T P6 = (-2,-2,0). (2.10) 



8-point function: 



The momenta are 



yo=\vm\, r=J(l-y*)(l-y*). (2.11) 



2n Pl = (-1,0,-1), 2n P2 = (1,0,-1), 2vrp 3 = (-1,1,0), 

2 m = (1,1,0), 2n P5 = (-1,0,1), 2vrp 6 = (1,0,1), 

2vrp 7 = (-1,-1,0), 27rp 8 = (1,-1,0). (2.12) 

These solutions do not reproduce the BDS formula. The most significant discrepancy 
appears in their infrared singularities. In fact, these solutions have the same infrared 
divergences as the 4-point amplitude. These solutions have also the delta function source 
term in the equation of motion. It has been discussed in [9] that these are not the minimal 
surface and correspond to other disconnected diagrams in field theory. 

The main purpose of this paper is to explore numerically the minimal surface corre- 
sponding to the same boundary conditions as the above higher-point solutions. In this 
paper we will study the square worldsheet for simplicity. We will make a comment on a 
generalization to the light-like hexagon Wilson loop in Sect. 4. 
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2.2 Discretized Euler-Lagrange Equations 

We firstly introduce the square lattice with spacing h — where M is a positive integer. 
At each site (—1 + hi, — 1 + hj) (i, j — 0, • • • , M), we assign the variables 

y [i,j]=y (-l + hi,-l + hj), r[i,j]=r(-l+hi,-l+hj). (2.13) 

We use the central difference method to discretize the Euler-Lagrange equations. Namely, 
for a function f(yi, y 2 ) we adopt the following rules to obtain the difference equations from 
the differential equations: 

f[i + lj]-f[i-lj] 



dif(yi,y 2 ) 
d 2 f(yi,y 2 ) 
did 2 f(yi,y 2 ) 



2h 

f[j,j + l]-f[j,j-l] 
2h 

f[i + 1, j + 1] - f[i + 1, j - 1] - f[i - 1, j + 1] + f[i -hj- 1] 

Ah 2 

f{ l + l,j]-2fii,j] + f{i-l,j] 
h 2 



alf{yuV2) _ /Miil-VMHffij-H (214) 

We need to specify the boundary conditions to solve the equations. For example, the 
boundary conditions (12.51) for the 4-point solutions lead to 

2/oM] = yo(-l+hi,-l), y [i, M) = y (-l + hi, 1), 
yo[QJ] = y (-h-l + hj), y [M,j] = y (l,-l + hj), 
r[i,0] = r[i,M}=r[0,j] = r[M,j]=0. (2.15) 

Then we obtain 2 x (M — l) 2 nonlinear simultaneous equations for yo[i,j] and j]. 
We will use Newton's method to find a numerical solution. In this method the initial 
solution is important to get good numerical results iteratively. The approximate solution 
will soon converge if the initial numerical data is appropriate. Otherwise, the numerical 
solution does not converge and the surface would not be smooth. We take as the initial 
condition the Alday-Maldacena solution for 4-point function or cut and glue solutions for 
6 and 8-point functions. But some trials show that the cut and glue solution does not 
lead to convergence when the lattice size becomes large. So we begin with the M — 10 
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lattice with these initial conditions and proceed step by step to larger size of lattices, up to 
M = 52cQ, with the output of a previous smaller lattice calculation being the input for the 
next larger lattice, by linearly interpolating the output. For each lattice size, the Newton's 
method is repeatedly applied until the non-linear simultaneous equation is satisfied up to 
O(l0~ 16 ) and the area of the obtained surface does not change up to O(10~ 6 ) even when 
we proceed to the next step of the Newton's method. The area is approximately evaluated 
as S = ^L[i,j]/i 2 , where L[i,j] and h are the discretized Lagrangian at a lattice point 
and the lattice spacing, respectively. 
We plot numerical solutions at M = 50 in Fig. HJ where we can see the 4-point solution 
apparently agrees with the exact solution. This agreement with the exact solution is 
quantitatively checked by comparing the areas between numerical and exact solutions in 
the next section. In Figs. [2} [3] and HI we plot numerical solutions of the 6-point and 8-point 
amplitudes with the same boundary conditions as the cut and glue solutions. We find 
that the numerical solutions are different from the cut and glue solutions and new cusps 
seem to appear which are absent in the cut and glue solutions \9\. In the next section we 
will evaluate the area of the surface numerically and examine the infrared behavior of the 
regularized area. 




(a) (b) 
Figure 1: Numerical solution of the 4-point amplitude at M = 50. 



1 A numerical computation was carried out by Mathcmatica and C programs on a personal computer 
with 8GB memory, which restricts the lattice size to M = 520. 
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(a) (b) 
Figure 2: Numerical solution of 6-point amplitude solution 1 (a) y (b) r at M = 50. 




0.8 - 

0.6 - 

0.4 - 

0.2 - 




0.8 - 

0.6 - 

0.4 - 

0.2 - 



1.2 r 
1 - 




M=50 
eq. (2.11) 



(a) 



(b) 



Figure 4: Numerical solution of 8-point amplitude (a) yo (b) r at M = 50. 

3 The cut-off regularization and the minimal area 

In order to study the gluon scattering amplitude, we need to evaluate the area of the 
minimal surface in AdS spacetime. Since we have obtained numerically the solution 
of the surface from the discretized Euler-Lagrange equations, we can evaluate it easily. 
But for a light-like Wilson loop, the value of the area diverges due to the cusp [HI [T]. 
We often use the dimensional regularization scheme to control this divergence, which is 
convenient to compare the area with the field theory result. But in order to calculate 
the area numerically in the present discretization procedure, it is better to use cut-off 
regularization. In a recent paper [13] . a cut-off r c in the radial coordinate is introduced 
to regularize the action, which characterizes the cusp divergences in a simple way. In this 
section, we will compare our numerical results with the cut-off formula and see that our 
results numerically agree with the cut-off regularized formula of the n-point amplitudes. 

3.1 The worldsheet cut-off regularization 

Before going to the radial cut-off formula, we will consider a simple worldsheet cut-off 
regularization, where we restrict the integration region of (y±, 2/2) i n the action (12.21) to 
[— 1 + 5, 1 — 5} x [— 1 + 5, 1 — 5} for small 5. Substituting the 4-point solution (12.61) into 
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the Lagrangian L, we get 



Here we omit the overall factor -f- 



L 




defined by 



S[S] = / dy 1 dy 2 

J -1+8 




(l-2/?)(l-2/2 2 )' 



1 



(3.2) 



This integral is evaluated easily and we obtain 



S[*]=(bg^) 



(3.3) 



The infrared singularity around the cusp appears as the (log5) 2 divergence, which corre- 
sponds to the double pole \ in the dimensional regularization scheme. 

We have seen that S[8] for cut and glue solutions of 6 and 8-points amplitudes is 
the same as the 4-point amplitude|9j. In Fig. we plot the 5 — S[5] graph for 4, 6, 
and 8-point functions. From the graphs we see that there exist differences among the 
regularized actions, especially near the boundaries (i.e. small 5). This fact indicates that 
our numerical 6 and 8-point solutions have different divergent properties from cut and 
glue solutions. We see differences between numerical and exact results for the 4-point 
amplitude at small 5 region. These are due to the fact that the numerical results are 
quite sensitive with the errors which come from the discretization of r near the boundary. 
However, the discrepancy at a given 5 is suppressed when the lattice size increases, as we 
will see later. 

Since we do not know exact functional form of new numerical solutions, the ^-dependence 
of the numerical S[S] for higher-point amplitudes is not clear. Roughly we can identify 
log 5 as -. It would not be a bad approximation to use a fitting function: 



for S[5}. Table [T] is a list of the value of a, b and c for each data fitted in the region 
0.04 < 5 < 1.14. The leading singularities are different for 4, 6 and 8 point amplitudes. 
But for solution 1 and solution 2 of 6-point functions they agree: a 4pt < aeptsoii ~ a 6ptsoi2 < 



a(log<5) 2 + 6(log5) + c, 



(3.4) 
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Figure 5: S - S[5] graphs at M = 300. 





a 


b 


c 


4-point 


1.12352 


-1.16767 


0.101811 


6-point sol. 1 


1.40647 


-0.648713 


0.214939 


6-point sol. 2 


1.41405 


-0.58385 


0.241626 


8-point 


1.67364 


-0.26576 


0.280678 



Table 1: Fitting parameters of M — 300 data using ( 13. 41) . 



10 



a8pt- This is consistent with the expected infrared leading singularity behaviour const, x 
n(log<5) 2 for the n-point amplitude [2]. 

It is a hard problem to estimate rigorously numerical errors in our approximation. 
But in the case of the 4-point amplitude, for which the exact solution has been obtained, 
we can see several positive features which support the validity of our numerical approach. 
Fig. [6] shows the behaviour of the numerical solution for r on the y% = plane near 
the boundary. We can see that the numerical solution approaches the exact one as M 
increases. Although we omit further details, it can be also shown that the solution is 
not affected by smooth and small modification of the initial condition. In addition, the 
area of the surface approaches to the exact result as M becomes larger and larger. In 
Table [2J and Fig. [TJ we exhibit the behaviours of the cut-off regularized area for 5 = 0.1 
and 8 = 0.2, which shows good convergence to the analytical results, with the error being 
about 0.2% for the 5 = 0.2, M = 520 case. 

0.2 > 
0.18 " 
0.16 
0.14 
0.12 

0.1 

0.08 

0.06 

0.04 

0.02 


0. 



i 1 — 

Exact sol. - 
M=100 - 
M=300 -— x- 
M=520 b 




Figure 6: Near boundary behavior of r on yi = plane for several M's (4-point solution). 

This artificial worldsheet cut-off is not invariant under conformal transformation of 
the worldsheet. This regularization heavily depends on the static gauge. In order to 
regularize the action in a conformal invariant way, it is better to introduce the radial 
cut-off as in [T3] . 
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Figure 7: Regularized area at 5 — 0.1 and 5 = 0.2 for several M's (4-point solution). 
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M 


5*4 


(5*4 — Se)/ S e 


100 


8.97215 


0.03488390 




100 


4.90702 


0.01641050 


200 


8.80018 


0.01504770 




200 


4.86233 


0.00715303 


300 


8.75113 


0.00939031 




300 


4.84934 


0.00446242 


520 


8.71324 


0.00501923 




520 


4.83927 


0.00237603 




(a) S 


= 0.1 




(b)5 


= 0.2 



Table 2: Convergent behavior of cut-off regularized area for the 4-point solution. S4 and 
are numerical and exact results, respectively. 
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3.2 The radial cut-off regularization 

In the radial cut-off regularization scheme we introduce a cut-off r c in the radial direction [T3 
The regularized area is surrounded by the cut-off curve C : r c = r(y 1 ,y 2 ). For example, 
the curve of the 4-point amplitude with s = t is 

r2 = (l -!/?)(! (3-5) 
The regularized action is the area of the region S whose boundary is C: 

^'/.^ rang - (3 ' 6) 

For fixed yi, the integration region of yi is — y\ < y<i < y%, where y%(> 0) is given by 

(Z/ 2 C ) 2 = 1 " j-^- (3-7) 
After the integration over y%, we get 



SH = rf * log f ^Zg + V^glg | . ,3.8) 

Since we are interested in the small r c limit, we expand the integrand f(yi,r c ) in r c : 
f ft/i r ) - 21ogr e -log(4-4y?) r e 2 

The integral over y\ leads to 

5[r c ] = ~ (log ^ 2 - 2(log2) 2 - ^ + 0(r c 2 ). (3.10) 

The general (s, i) solution has been found in [I], which can be obtained by scale and 
boost transformation of the s = t solution: 

r /_ ar , _ aVT+¥yo , _ ayj , , 

1 + 6y 1 + &i/o 1 + %o 

where a is a parameter for the scale transformation and b the boost parameter. The 

Mandelstam variables s and t are given by 

W« = -(i^. ^) 2t = -(T^- < 3 - 12 > 
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Figure 8: (a)The cut-off curves with — = 0.3 and 6 = 0, 0.5, 0.9 (b) The cut-off curve 
with — = 0.3 from the numerical data (4-point solution M = 300). 

The cut-off curve (I3.5P (see Fig. [8] (a)) is now replaced by 

^(l-tfU-tfj—L-p. (3.13) 

We can put a = 1 by rescaling r c — > r c a. For fixed y±, y 2 takes the value in the range 
yf <V2< y 2 + , where 

, (C ± _ -brfyi ± V(l ~ yl){l ~ y\ - r 2 c + b 2 r 2 y 2 ) 

V2 ~ l-y 2 + b 2 r 2 c y 2 ' ^ i4j 

The integral over y 2 yields 



$[>:}= I" i-^-dyif(yi,r c ), (3.15) 

l-b 2 r'£ 



where 

1 1, fl + y° + l-y c 2 



S^—^-^-^y ,3, 6) 

Expanding in r c we get 

ffarc) = -T^^W^)^ ^ (3 ' 1T) 
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After the integral over t/i, we get 



S[r c ] = - log 2 c — + - log 2 — V - - log R - 2(log 2Y + Oiri). 

(3.18) 

This formula of the 4-point amplitude was obtained in [T3] using the conformal gauge. 

The above analysis for the regularized action of the 4-point amplitude is generalized 
to examine the infrared singular part of the n-point amplitude, which was done in the 
dimensional regularizationJB]. The infrared singularity of the n-point amplitude is char- 
acterized by the cusp made of two external gluon lines with momenta pi and p%+i. It is 
convenient to use the light-cone coordinates y = y° ± y 1 and regard y 2 and r as the 
functions of y ± . The Nambu-Goto action in this gauge is 

hi 



S = — / dy_dy + —y/l - Ad-y 2 d + y 2 - 4<9_r<9 + r - A(d^y 2 d+r - d-rd + y 2 ) 2 . (3.19) 



The momenta of the external gluon lines in the (y-,y + ,y 2 ) coordinates are parametrized 
as 27ipi = Z\(0, 1, 0) and 2np i+ i = z 2 (a, 1, \fa). The Mandelstam variable is (27r) 2 Sj 5 j + i = 
Z\z 2 (—a) = —z\Z 2 ol . Then the solution of the equation of motion, which approaches to 
the cusp solution in the limit a — > 0, is 



r(y-,y+) = v^y y- (y+ - ^y^j , y 2 (y-, y+) = -^V- ( 3 - 20 ) 
Parameterizing the solution as 

y _ = az 2 Y_, y + = Zl Y+ + z 2 Y„, (3.21) 

we find 



r(Yl,r+) = v / 2 v /y_F + y / -(27r) 2 SM+1 , (3.22) 
and the action near the cusp is 

<W = £f£^. (-3) 



4tt J 2Y_Y 

The action is divergent at the cusp Y_ = Y + = 0. Introducing the radial cut-off r c by 



v^ v /-(27r)2 SM+1 ^K3^, (3.24) 
15 



the regularized action rescaled by the factor of j- becomes 

1 f l 



i Z" 1 „, f 1 i 



2 Jr§ ^7 rl 2YY 

- K l0g ^fcr) 2 - (3 ' 25) 

Here A = y/2y/-(2ir) 2 s iii+ i. 

Then the n-point amplitude is expected to have the structure 

= T kg + F„(p!, • • • , Pn ) + 0(r 2 ), 



(3.26) 



where s^+i = —{pi + Pi+i) 2 and p n+ i = p\. F n (pi, • ■ ■ ,p n ) is a finite remainder part. 
The 4-point function with the momenta (12.41) with s = t, we obtain (27r) 2 si.2 = —8. The 
action is of the form 2 

^N = y(log^) +c , (3.27) 

up to 0(r 2 ) terms, where a = 1 and Co is a constant. 

For the 6-point function (solution 1 and 2) and the 8-point solution, the amplitudes 
in the radial cut-off scheme are 

Sl[r c ] = ^ [4 Tlog 2 + ("log 2 + ("log [|) 2 ) + ci, (3.28) 



St[r c ] = ~a 2 (log^) +c 2 , (3.29) 



2\ 2 



Ss[r c ] = log -M + log-r + c 3, (3-30) 




up to 0(r 2 ) terms, where Cb\ — CL2 — O3 — 1 and ci, C2 and C3 are constants. 

We approximate the area by summing the value of discretized Lagrangian at the lattice 
points at which the value of r[i,j] is larger than the radial cut-off r c . In Fig. [9] we 
show the fitting data of the r c regularization. The numerical results are fitted to the trial 
functions with the least-square method in the range 0.15 < r c < 3.0. The lower limit of 
the fitting range cannot be so small because of the rapid growth of r near the boundary 
(see Fig. E]). We reasonably chosen the lower limit as r c = 0.15 so that the lattice points 
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around the boundary which satisfy r[i,j] > r c are several points inner from the boundary. 
The approximation error for the area is roughly estimated as the difference between the 
area estimated above and the one for the nearest outer lattice points. Fig. [10] is the result 
with such errors for the 4-point solutions. We fit the results S[r c ] for each solution with 
the weight of errors estimated in this way by using 'fit' command of GNUPLOT. Fitting 
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Figure 9: r c - S[r c ] graph at M = 520. 

these data in the region 0.15 < r c < 0.3 by the functions (I3T27D . fl3T2^|) and (137501) 
respectively, we find that 

a = 1.01694, c = -3.5912, (3.31) 

a x = 1.03715, ci = -3.68904, (3.32) 

a 2 = 1.0376, c 2 = -3.70123, (3.33) 

a 3 = 1.04973, c 3 = -3.70072. (3.34) 

We see that our numerical results are in agreement with the radial cut-off regularization 
formula within a few per cent. The constant Co of the 4-point amplitude differs from the 
cut-off formula (13. 10p . Therefore we need to take into account for the finite r 2 c corrections 
to the area in order to get the finite remainder part. This subject is beyond the scope of 
the present paper and is left to a future problem. 
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Figure 10: r c - S[r c ] graph at M — 520 with error-bars. 



4 Conclusions and Discussion 

In this paper we have investigated the discretized Euler equations of the minimal surface 
in AdS spacetime numerically. We examined the minimal surface corresponding to the 4, 
6 and 8-point amplitudes, in which the n > 6 solutions have the same boundary conditions 
as in [9|. These boundary conditions are simple because the worldsheet in the static gauge 
is square, which are easy to write a program to solve the system of nonlinear simultaneous 
equations. Since the area is divergent around the cusps, we have introduced the worldsheet 
cut-off and radial cut-off regularizations. In the worldsheet cut-off regularization, the 4- 
point amplitude agrees numerically with the exact solution within 0.2% for the 5 = 0.2, 
M = 520 case. For the radial cut-off regularized amplitudes, we have checked that the 
infrared singularity part of the area numerically agrees with the analytical result near the 
cusp within 5% for 520 x 520 lattice points. 

Using conformal transformation of the solution, we are able to study non-trivial mo- 
mentum dependence of the amplitude including not only the infrared singularity but also 
the finite part. It is an interesting problem to extract numerically the finite remainder 
part of the gluon amplitudes from the minimal surface. 

One can also generalize the present square worldsheet to polygons. Here we will show 
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an preliminary example of the minimal surface solution with the hexagonal boundary 
(Fig. II ip . whose momenta are 

2rr Pl = (V2,l,-1), 2irp 2 = (-V2, -1, -1), 2np 3 = (1, -1, 0), 

2ttp4 = (-V2,-l,l), 2 m = ( v / 2,l,l), 2 m = (-1,1,0). (4.1) 



Here we can see six cusps. The evaluation of the regularized action with polygon light- 
like boundary condition gives another nontrivial test of the gluon amplitudes/ Wilson loop 
duality at strong coupling. 




Figure 11: Numerical solution of 6-point function with hexagon boundary (116 lattice 
points) (a) y (1/1,3/2), (b) r(y 1 ,y 2 ). 

Another interesting application of this formalism is the gluon scattering in non-AdS 
geometry, which is rather difficult to obtain the analytic solution of the Euler-Lagrange 
equations jTTJ. A numerical approach will help us to understand their properties at strong 
coupling. A detailed study of these problems will be discussed elsewhere. 
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